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Abstract. The interpretation of the difference observed between the 
radial distribution of globular clusters and that of halo-bulge stars in 
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elliptical galaxies is discussed in terms of evolution of their globular clus- 
ter systems. I present a short summary of the evidence that dynamical 
evolution of globular cluster systems is not only able to explain the flat- 
tening of their distribution toward the galactic center, but also it may 
have played an important role on the primordial activity of the parent 
galaxy. 



1. Introduction 

o 

, It is well known that the haloes and bulges of galaxies constitute 'collision- 

Qh' less' systems, for its 2-body relaxation times are mich longer than their age. 

q Globular clusters (GC), too, are halo objects but they are characterized by a 

mass 10 5 — 10 6 heavier than the stars of the background where they orbit: this 
implies they do not constitute a collisionless system. 

The collective effect of the fluctuating gravitational field produced by light 
stars along the cluster orbital motion reflects in the well known dynamical fric- 
tion braking, that has been shown by various authors to be particularly impor- 
tant in triaxial galaxies, where the lack of symmetries in the potential excludes 
the conservation of any component of the angular momentum. This frictional 
braking acts on a time scale inversely proportional to the cluster mass and in 
triaxial galaxies it can be an order of magnitude shorter than in symmetric 
galaxies. Moreover, for sufficiently massive clusters, the orbital energy is lost in 
less than 1 Gyr. In addition to the dynamical friction collective effect, on glob- 
ular clusters acts also a local process: the tidal interaction with the (compact) 
galactic nucleus. The 2 evolutionary processes compete, in some sense, because 
dynamical friction tends to concentrate the GCS distribution around the galac- 
tic center, while a massive galactic nucleus (a super massive black hole?) may 
destroy a globular cluster which approaches too much to it, and this corresponds 
to an erosion of the GC inner galactic population. 

The balance of the 2 mechanisms may lead to the explanation of the difference 
observed between the GCS and the halo-bulge stellar component radial distri- 
butions in some galaxies: the globular cluster systems (GCSs) are sistematically 
less concentrated to the parent galaxy center than bulges. Another important 
consequence of the evolution of the GCS in a galaxy is the quantity of mass 
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that, in form of stars tidally stripped by the nucleus, can feed the massive cen- 
tral black hole, releasing gravitational energy that can be radiated in form of 
electromagnetic waves, sustaining a violent nuclear activity. In this brief report 
we will describe and discuss some relevant aspects of this fascinating scenario, 
including that concerned with the stellar population of the inner kiloparsec of a 
galaxy 

2. The globular cluster distribution in a galaxy 

The comparative discussion of the radial distributions of the stellar component 
and of globular clusters in elliptical galaxies has started not much time ago, due 
to the difficulty of determining in a reliable way the globular cluster distribution 
of the central galactic regions, even in nearby galaxies. Grillmair et al. (1986) 
were the first to show how the stellar bulge in M 49 is more concentrated to the 
center of the galaxy than the GCS one. Much other evidence of such phenomenon 
was later given (see, for instance, Lauer & Kormendy 1986; Harris 1986; Harris 
et al. 1991), and now we can say that no case has been found yet where the 
GCS distribution in a galaxy is more peaked to its center than the surrounding 
halo-bulge stellar population. 

The observational evidence does not have a certain interpretation. 
Some claim that the reason of the different profiles is due to a difference in the 
formation epoch of the 2 halo components: globular clusters formed earlier, when 
the protogalaxy was more expanded, than the collisionless stellar component. 
Some others, starting from the reasonable hypothesis that the various com- 
ponents of the halo population formed almost simultaneously, think that the 
present difference is due to an evolutionary process. No clear evolutive processes 
can affect the stellar component as a whole, while the globular cluster distri- 
bution in a galaxy can suffer for both the collective effect of dynamical friction 
induced by stars and for the tidal, quasi impulsive, interaction with a compact 
nucleus sited in the galactic center. Shortly, dynamical friction is more effective 
on heavier clusters that are, usually, more tightly bound and so resist better to 
tidal shocks. As a consequence, the radial distribution of globular clusters in a 
galaxy should be mass (and density) dependent, provided the time scales of the 
mentioned effects are short enough in comparison with the age of the system. 
We do not go here into details but just mention the work by Pesce, Capuzzo- 
Dolcetta &: Vietri (1992) and Capuzzo-Dolcetta (1993) who quantitatively showed 
how relevant are the frictional and tidal effects on the evolution of a GCS, partic- 
ularly in a triaxial galaxy. Probably, the most relevant outcome of these works 
is that the time scales of the 2 evolutionary processes can be (in dependence on 
mass and density of clusters, indeed) short enough to be surely relevant for the 
evolution of GCSs from their initial state. 

A straightforward application of the results of Pesce et al. (1992) and 
Capuzzo-Dolcetta (1993) lead Capuzzo-Dolcetta & Tesseri (1997) to a succesful 
comparison of the evolved theoretical GCS spatial distributions (at varying some 
relevant parameters of the model) with the general characteristics of the observed 
distributions like, for example, the observed larger core radii of the GCS radial 
distributions with respect to the bulge one. So, the evolution of the GCS in 
an elliptical galaxy seems a confirmed hypothesis, and, other explaining the 
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mentioned difference in radial profiles, it may have played an important role in 
the primeval lives of galaxies, supplying a significant amount of matter to the 
central galactic regions, often sites of violent activity (AGNs). 

2.1. The comparison between the GCS and stellar bulge radial dis- 
tributions 

In addition to our galaxy, M 31 and M 87 (discussed in Capuzzo-Dolcetta Sz 
Vignola 1997) the comparison between the stellar bulge and globular cluster 
radial distributions has been extended in details to a set of 14 elliptical galaxies 
observed with the WFPC of the HST by Forbes et al. (1996), Forbes et al. 
(1998a,b). The results, presented in Capuzzo-Dolcetta & Tesseri (1999) and 



Table 1. The presently observed number of clusters (N), its initial 
value (Ni), its fractional variation (A), the mass lost in form of dis- 
appeared globulars (Mi). Data are from Capuzzo-Dolcetta & Vignola 
(1997), Capuzzo-Dolcetta & Tesseri (1999) and Capuzzo-Dolcetta & 
Donnarumma (2000). 



Galaxy 


N 


Ni 


A 


M t (M ) 


Milky way 


155 


211 


0.26 


1.80 x 10'' 


M 31 


283 


368 


0.23 


2.30 x 10 7 


M 87 


4456 


8021 


0.44 


2.33 x 10 9 


NGC 1379 


132 


512 


0.74 


1.50 x 10 s 


NGC 1399 


5168 


9680 


0.63 


1.44 x 10 s 


NGC 1404 


508 


1061 


0.53 


1.75 x 10 s 


NGC 1427 


248 


487 


0.49 


8.86 xlO 7 


NGC 1439 


130 


141: 


0.08: 


4.79 xlO 6 : 


NGC 1700 


25 


39: 


0.36: 


3.66 xlO 6 : 


NGC 4365 


517 


849 


0.39 


7.48 x 10 7 


NGC 4494 


200 


297 


0.33 


2.98 xlO 7 


NGC 4589 


241 


371 


0.35 


7.58 xlO 7 


NGC 5322 


175 


266 


0.34 


6.51 xlO 7 


NGC 5813 


382 


596 


0.36 


1.03 xlO 8 


NGC 5982 


135 


260 


0.48 


8.86 xlO 7 


NGC 7626 


215 


365 


0.41 


3.59 xlO 8 


IC 1459 


271 


516 


0.47 


1.57 xlO 8 



Capuzzo-Dolcetta & Donnarumma (2000), show in all the cases a significantly 
flatter distribution of the GCS. 

Assuming that initially the stellar bulge and the GCS had the same profile it 
is easy to evaluate the number of globular clusters 'lost ' in a galaxy as the 
difference of the surface integrals of the 2 profiles. A proper estimate of the 
mean value of the globular cluster initial mass brings to a value of the quantity 
of mass lost to the central galactic region. This was done in the mentioned 
papers; as a summary of the results, in Table 1 we give the values relative to 
the set of 17 galaxies studied so far. It is really remarkable the high percentage, 
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41%, of the initial GCS mass that, according to our estimates, is fallen, in the 
average, in the inner galactic zones. The quantity of mass lost to the central 
regions, averaged over the set of 17 galaxies, is 2.19 x 10 8 M , that should have 
been relevant as a source of fuel for a possible central compact object. 

3. Nucleus feeding by centrally merged globular clusters 

As shown in Capuzzo-Dolcetta (1993) the possible causes of the significant dif- 
ference observed between the stellar bulge and the GCS distribution in a galaxy 
are dynamical friction and tidal disruption, this latter induced by a compact 
galactic nucleus, if present. In that paper a preliminary quantitative discus- 
sion of the modes of formation of a supercluster composed by orbitally decayed 
clusters merged at the galactic center, as well as of debris of tidally shattered 
ones has been discussed and demonstrated to occur in times acceptably short 
to give birth to a supermassive object which could release stars to the central 
black hole in a way to make it a powerful engine of electro-magnetic energy 
extracted by the gravitational field. A subsequent, still in preparation, paper by 
Capuzzo-Dolcetta (2001) presents a detailed quantitative model of the evolution 
of a GCS in a triaxial galaxy leading to a detailed estimate of the time growth of 
the central compact seed as well as of the luminosity evolution of the resulting 
AGN. Of course, the model depends on some parameters like, mainly, the initial 
velocity dispersion of globular clusters, their IMF, their orbital distribution, etc.; 
these parameters can just partially be fixed a priori. By the way, it is relevant 
noting that a limitation of the possible values of these parameters is induced 
by the request that the models give, at the present time, acceptable numbers 
of observable clusters with the proper observed radial distribution and velocity 
dispersions, in order to be fully self-consistent. 

We do not go here into details of the models, but rather want just to give a 
qualitative insight to the phenomenology, starting from first principles, that is 
useful to convince that the overall scenario is realistic at all. 
Suppose that around a central galactic black hole (b.h.) of mass there is 
an almost steady distribution of stars with typical mass and radius m* and i?*, 
and having space density and average velocity < >, then the rate of mass 
accretion onto the b.h. is simply m = p*a < v* >. In the mass accretion rate a 
is the 'capture ' cross-section given by 

( Gm bh \ 
i+ r^)=" 3 ( 1+a ^ £ )' 

where ri n f is the 'influence ' radius of the b.h. and is the 'destruction ' radius 
(below this distance from the b.h. the star is tidally destroyed or swallowed) 
given by 

r d = max(r t ,rs), 

where rj is the 'tidal ' radius and r$ the Schwarzschild's radius of the b.h., 
respectively. Because the ratio 

2 I 

r s _ 2AGm'i § 
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Figure 1. panel a: Bi-logarithmic plot of the curve (solid line) 
p* < > _1 (in units M pc _3 /(kms -1 )) that corresponds to rh = 1 
M yr _1 in function of mbh- The dashed line is the 'collisional ' curve 
(see text). 

panel b: The correlation between the (logarithmic) GCS mass lost and 
the central galactic black hole mass for 7 galaxies for which these data 
are available. Masses are in solar masses. The straight solid line is the 
least square fit to the data. 



is > 1 just for massive b.h. (m^ > m^crit — 1-62 x 10 8 M Q assuming solar 
values for the cluster star mass and radius) it is clear that the process of b.h. 
feeding is by mean of tidally disrupted stars until the b.h. mass has grown enough 
to swallow the surrounding stars whole when entering the event horizon through 
the Schwarzschild's radius. In the 2 cases we have (assuming < vl >=< > 2 ) 



m(M Q yr- 1 ) = { 8 " 2 X 10 l^th™* ^ ''P* R * < v * > \ r d = r t , 
V ' \ 2.8 x 10- 21 ml hP * < >- x r d = r s , 

respectively (m*, m^, p* and i?* are in solar units and < > in km s" 1 ). Let 
us note that, of course, the efficiency of conversion into radiation is higher in 
the case when stars are tidally destroyed before falling onto the compact object 
than when swallowed whole. 

On the basis of previous equations, it is possible to draw the plot given in Fig. 
la, where, in function of m^, the minimum values of < v* >~ x required to 
yield an accretion rate of at least 1 M Q yr _1 are given (1 M yr _1 is the typical 
value of mass accretion required to sustain an AGN) yielding a separation into 
a 'high' accretion and a 'low' accretion density zone. 

As an example, the known values of the central b.h. masses in our galaxy, 
M31 and M 87 (2.6 x 10 6 , 6.2 x 10 7 , 3.6 x 1O 9 M , respectively) imply a stellar 
density at least 4 x 10 7 M Q , 6.3 x 10 5 M , 400 M , respectively, to accrete 1 M 
yr _1 in the hypothesis of < t>* >= 10 km s _1 (which is, anyway, just a typical 
value of the mean stellar velocity in a globular cluster). It is remarkable that 
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these density values are not too high to be reached by a stellar system composed 
by merged globulars, at least in its central region (many galactic globulars have 
central densities greater than 10 4 M pc~ 3 , and few of them greater than 10 5 
M.q pc~ 3 ). It may be interesting noting that in particular situations of very high 
dense clusters surrounding the compact nucleus, the accretion can occur on the 
larger space scale given by the collisional radius 

_ Gm bh _ 1 m bh 

r C0ll 9 r. -K*J 

(w e * is the escape velocity from star surface) which determines the radius of 
the sphere centred on the b.h. where star-star physical collisions are disruptive 
(see e.g. Binney and Tremaine 1998). Actually, it turns out that in normal 
conditions r co u is the largest of the 3 radii (rj, rs, r co u), being 

' 2^ « 1, 



^coll Ft* 

being rs* is the star Schwarzschild radius, and 



n 



= 2 4/3 f^f 3 ^ 



r co ll \m bh 



whenever m^/m* > 4. 

A larger corresponds, of course, to a larger a and so it would allow a huge 
mass accretion onto the b.h., as given by 



rh coll (M Q yr x ) = 3.27 x 10 16 m^m, >*#* < v* > 1 • 

For the sake of comparison, the line corresponding to this physically collisional 
regime is also shown in Fig. la. 

To conclude, it is a relevant point in favour of the GCS evolutionary picture 
the existence of a positive correlation (as shown in Figure lb) among the values 
of the (possible) mass lost in form of orbitally decayed clusters (Mj in Table 1) 
and those of the central b.h. in the 7 galaxies for which these data are available 
(the least square straight line fitting the whole data set has a correlation factor 
r 2 = 0.61 and standard error = 0.56). Actually, as remarked in Capuzzo™ 
Dolcetta & Tesseri (1999) a positive correlation is expected because larger black 
holes are grown during the evolution and they should induce a greater evolution 
of the GCS. 



4. Conclusions 

I have presented a summary of some results concerning the globular cluster dis- 
tribution in galaxies as compared with that of the underlying stellar component. 

It is known that globular clusters have a less peaked distribution; if this 
centrally flatter behaviour is interpreted in evolutive terms, it may imply a 
number of 'missing ' clusters that can be estimated in the assumption of an 
initial cluster distribution equal to that of stars of the halo-bulge. The found 
number of globular clusters likely lost to the inner galactic region is evaluated 
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to be significant in the galaxies studied so far. This may have had important 
consequences on the initial life of the galaxy, for it corresponds to a large quantity 
of mass in form of stars that have been buzzing around the center of the galaxy 
feeding a compact nucleus therein and allowing it to emit energy at the levels 
of AGNs. Another consequence of the globular cluster system evolution in a 
galaxy is the contribution to the inner field, in a way that has not been yet 
deeply analysed. 

Rather than going into details of the model describing the evolution of the 
coupled globular clusters-central black hole system in an elliptical galaxy (that 
can be found in other longer papers referred in this report), in this short paper 
I preferred to give simple evaluations of the densities required by the cluster 
around the galactic central black hole in order to be able to release it a sufficient 
quantity of mass to feed it as engine of an AGN. 
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